Unravelling the maternal evolutionary history of the African leopard (Panthera pardus pardus).

The African leopard (Panthera pardus pardus) has lost a significant proportion of its historical range, notably in north-western Africa and South Africa. Recent studies have explored the genetic diversity and population structure of African leopards across the continent. A notable genetic observation is the presence of two divergent mitochondrial lineages, PAR-I and PAR-II. Both lineages appeared to be distributed widely, with PAR-II frequently found in southern Africa. Until now, no study has attempted to date the emergence of either lineage, assess haplotype distribution, or explore their evolutionary histories in any detail. To investigate these underappreciated questions, we compiled the largest and most geographically representative leopard data set of the mitochondrial NADH-5 gene to date. We combined samples (n = 33) collected in an altitudinal transect across the Mpumalanga province of South Africa, where two populations of leopard are known to be in genetic contact, with previously published sequences of African leopard (n = 211). We estimate that the maternal PAR-I and PAR-II lineages diverged approximately 0.7051 (0.4477-0.9632) million years ago (Ma). Through spatial and demographic analyses, we show that while PAR-I underwent a mid-Pleistocene population expansion resulting in several closely related haplotypes with little geographic structure across much of its range, PAR-II remained at constant size and may even have declined slightly in the last 0.1 Ma. The higher genetic drift experienced within PAR-II drove a greater degree of structure with little haplotype sharing and unique haplotypes in central Africa, the Cape, KwaZulu-Natal and the South African Highveld. The phylogeographic structure of PAR-II, with its increasing frequency southward and its exclusive occurrence in south-eastern South Africa, suggests that this lineage may have been isolated in South Africa during the mid-Pleistocene. This hypothesis is supported by historical changes in paleoclimate that promoted intense aridification around the Limpopo Basin between 1.0-0.6 Ma, potentially reducing gene flow and promoting genetic drift. Interestingly, we ascertained that the two nuclear DNA populations identified by a previous study as East and West Mpumalanga correspond to PAR-I and PAR-II, respectively, and that they have come into secondary contact in the Lowveld region of South Africa. Our results suggest a subdivision of African leopard mtDNA into two clades, with one occurring almost exclusively in South Africa, and we identify the potential environmental drivers of this observed structure. We caution that our results are based on a single mtDNA locus, but it nevertheless provides a hypothesis that can be further tested with a dense sample of nuclear DNA data, preferably whole genomes. If our interpretation holds true, it would provide the first genetic explanation for the smaller observed size of leopards at the southernmost end of their range in Africa.

NONDOMESTIC FELID ABC BLOOD PHENOTYPING, GENOTYPING, AND CROSSMATCHING.

Based upon previous clinical experience with domestic cats (Felis catus), the ability to assess ABC blood types and blood (in-)compatibilities of nondomestic felids, and adequately consider and plan for blood transfusions, may be important. Although nondomestic felids appear to have an ABC blood group system similar to domestic cats, typing with point-of-care kits and by CMAH genotyping for domestic cats have not been reported. In this study, 162 blood samples from 18 different nondomestic felid species (cheetah [Acinonyx jubatus, n = 42], lion [Panthera leo, n = 33], tiger [Panthera tigris, n = 23], Canada lynx [Lynx canadensis, n = 11], snow leopard [Uncia uncia, n = 10], puma [Puma concolor, n = 7], clouded leopard [Neofelis nebulosa, n = 6], serval [Leptailurus serval, n = 5], jaguar [Panthera onca, n = 5], fishing cat [Prionailurus viverrinus, n = 4], Pallas cat [Felis manul, n = 3], bobcat [Lynx rufus, n = 3], ocelot [Leopardus pardalis, n = 3], black footed cat [Felis nigripes, n = 2], leopard [Panthera pardus, n = 2], African wildcat [Felis lybica, n = 1], caracal [Caracal caracal, n = 1], and sand cat [Felis margarita, n = 1]) were ABC blood typed by laboratory and point-of-care tests, genotyped for four known CMAH variants for type B and type C (AB) phenotypes, and crossmatched with one another and domestic type A cats. Traditional tube typing identified blood type A (n = 106), type B (n = 8), type C (n = 43), and no discernible ABC type (n = 4). Several discrepancies were found between point-of-care and traditional typing test results. None of the tested felids possessed the four CMAH variants responsible for type B and C (AB) in domestic cats. Crossmatch incompatibilities (≥2+ agglutination) were identified within and between nondomestic felid species and beyond ABC incompatibilities. Of 26 crossmatches performed between domestic cats and various nondomestic felids, only 7 (27%) were compatible. In conclusion, point-of-care typing kits and CMAH genotyping, successfully used in domestic cats, may not identify the correct ABC blood type in nondomestic felids. Prior crossmatching is recommended to increase the likelihood of compatible transfusions between any nondomestic felids.

NEW RECORDS OF HELMINTHS OF THE JAGUAR IN MEXICO, WITH AN UPDATED LIST OF SPECIES IN THE AMERICAS.

An inventory of parasites infecting the jaguar (Panthera onca) across its distribution range is relevant for the conservation of this threatened big cat. In this study, we report the occurrence of helminths in a jaguar from Mexico using morphological techniques (cleared and stained mounts and scanning electron microscopy) and partial sequences of the 28S ribosomal RNA (28S rRNA) gene and the cytochrome c oxidase 1 mitochondrial (COI) gene. We also provide an updated list of helminth species reported in jaguars in the Americas. Three helminth taxa are identified in the jaguar examined from Mexico: Toxocara cati, Physaloptera sp., and Taenia sp. The new 28S rRNA sequences of To. cati, Physaloptera sp., and Taenia sp. and the COI sequence of Taenia sp. corroborate the identity of the helminths isolated from this host. One hundred and twenty-nine records of helminths parasitizing jaguars from 49 studies up to May 2023 were identified in the Americas. In most of these studies (73.6%), helminths were identified using coproparasitological techniques. Sixteen helminths (7 nematodes, 5 cestodes, 3 acanthocephalans, and 1 trematode) were identified at the species level in free-ranging and captive jaguars. The study demonstrates the value of an integrative taxonomy approach to increase the accuracy of parasite identification in wildlife, especially when helminth specimens are scarce or poorly fixed.

Diversity of selected toll-like receptor genes in cheetahs (Acinonyx jubatus) and African leopards (Panthera pardus pardus).

The anthropogenic impact on wildlife is ever increasing. With shrinking habitats, wild populations are being pushed to co-exist in proximity to humans leading to an increased threat of infectious diseases. Therefore, understanding the immune system of a species is key to assess its resilience in a changing environment. The innate immune system (IIS) is the body's first line of defense against pathogens. High variability in IIS genes, like toll-like receptor (TLR) genes, appears to be associated with resistance to infectious diseases. However, few studies have investigated diversity in TLR genes in vulnerable species for conservation. Large predators are threatened globally including leopards and cheetahs, both listed as 'vulnerable' by IUCN. To examine IIS diversity in these sympatric species, we used next-generation-sequencing to compare selected TLR genes in African leopards and cheetahs. Despite differences, both species show some TLR haplotype similarity. Historic cheetahs from all subspecies exhibit greater genetic diversity than modern Southern African cheetahs. The diversity in investigated TLR genes is lower in modern Southern African cheetahs than in African leopards. Compared to historic cheetah data and other subspecies, a more recent population decline might explain the observed genetic impoverishment of TLR genes in modern Southern African cheetahs. However, this may not yet impact the health of this cheetah subspecies.

A cost-effective blood DNA methylation-based age estimation method in domestic cats, Tsushima leopard cats (Prionailurus bengalensis euptilurus) and Panthera species, using targeted bisulphite sequencing and machine learning models.

Individual age can be used to design more efficient and suitable management plans in both in situ and ex situ conservation programmes for targeted wildlife species. DNA methylation is a promising marker of epigenetic ageing that can accurately estimate age from small amounts of biological material, which can be collected in a minimally invasive manner. In this study, we sequenced five targeted genetic regions and used 8-23 selected CpG sites to build age estimation models using machine learning methods at only about $3-7 per sample. Blood samples of seven Felidae species were used, ranging from small to big, and domestic to endangered species: domestic cats (Felis catus, 139 samples), Tsushima leopard cats (Prionailurus bengalensis euptilurus, 84 samples) and five Panthera species (96 samples). The models achieved satisfactory accuracy, with the mean absolute error of the most accurate models recorded at 1.966, 1.348 and 1.552 years in domestic cats, Tsushima leopard cats and Panthera spp. respectively. We developed the models in domestic cats and Tsushima leopard cats, which were applicable to individuals regardless of health conditions; therefore, these models are applicable to samples collected from individuals with diverse characteristics, which is often the case in conservation. We also showed the possibility of developing universal age estimation models for the five Panthera spp. using only two of the five genetic regions. We do not recommend building a common age estimation model for all the target species using our markers, because of the degraded performance of models that included all species.

A unique single nucleotide polymorphism in Agouti Signalling Protein (ASIP) gene changes coat colour of Sri Lankan leopard (Panthera pardus kotiya) to dark black.

The Sri Lankan leopard (Panthera pardus kotiya) is an endangered subspecies restricted to isolated and fragmented populations in Sri Lanka. Among them, melanistic leopards have been recorded on a few occasions. Literature suggests the evolution of melanism several times in the Felidae family, with three species having distinct mutations. Nevertheless, the mutations or other variations in the remaining species, including Sri Lankan melanistic leopard, are unknown. We used reference-based assembled nuclear genomes of Sri Lankan wild type and melanistic leopards and de novo assembled mitogenomes of the same to investigate the genetic basis, adaptive significance, and evolutionary history of the Sri Lankan melanistic leopard. Interestingly, we identified a single nucleotide polymorphism in exon-4 Sri Lankan melanistic leopard, which may completely ablate Agouti Signalling Protein (ASIP) function. The wild type leopards in Sri Lanka did not carry this mutation, suggesting the cause for the occurrence of melanistic leopords in the population. Comparative analysis of existing genomic data in the literature suggests it as a P. p. kotiya specific mutation and a novel mutation in the ASIP-gene of the Felidae family, contributing to naturally occurring colour polymorphism. Our data suggested the coalescence time of Sri Lankan leopards at ~0.5 million years, sisters to the Panthera pardus lineage. The genetic diversity was low in Sri Lankan leopards. Further, the P. p. kotiya melanistic leopard is a different morphotype of the P. p. kotiya wildtype leopard resulting from the mutation in the ASIP-gene. The ability of black leopards to camouflage, along with the likelihood of recurrence and transfer to future generations, suggests that this rare mutation could be environment-adaptable.

A new multiplex qPCR assay to detect and differentiate big cat species in the illegal wildlife trade.

All species of big cats, including tigers, cheetahs, leopards, lions, snow leopards, and jaguars, are protected under the Convention on the International Trade in Endangered Species (CITES). This is due in large part to population declines resulting from anthropogenic factors, especially poaching and the unregulated and illegal trade in pelts, bones, teeth and other products that are derived from these iconic species. To enhance and scale up monitoring for big cat products in this trade, we created a rapid multiplex qPCR test that can identify and differentiate DNA from tiger (Panthera tigris), cheetah (Acinonyx jubatus), leopard (Panthera pardus), lion (Panthera leo), snow leopard (Panthera uncia), and jaguar (Panthera onca) in wildlife products using melt curve analysis to identify each species by its unique melt peak temperature. Our results showed high PCR efficiency (> 90%), sensitivity (detection limit of 5 copies of DNA per PCR reaction) and specificity (no cross amplification between each of the 6 big cat species). When paired with a rapid (< 1 h) DNA extraction protocol that amplifies DNA from bone, teeth, and preserved skin, total test time is less than three hours. This test can be used as a screening method to improve our understanding of the scale and scope of the illegal trade in big cats and aid in the enforcement of international regulations that govern the trade in wildlife and wildlife products, both ultimately benefiting the conservation of these species worldwide.

Numt Parser: Automated identification and removal of nuclear mitochondrial pseudogenes (numts) for accurate mitochondrial genome reconstruction in Panthera.

Nuclear mitochondrial pseudogenes (numts) may hinder the reconstruction of mtDNA genomes and affect the reliability of mtDNA datasets for phylogenetic and population genetic comparisons. Here, we present the program Numt Parser, which allows for the identification of DNA sequences that likely originate from numt pseudogene DNA. Sequencing reads are classified as originating from either numt or true cytoplasmic mitochondrial (cymt) DNA by direct comparison against cymt and numt reference sequences. Classified reads can then be parsed into cymt or numt datasets. We tested this program using whole genome shotgun-sequenced data from 2 ancient Cape lions (Panthera leo), because mtDNA is often the marker of choice for ancient DNA studies and the genus Panthera is known to have numt pseudogenes. Numt Parser decreased sequence disagreements that were likely due to numt pseudogene contamination and equalized read coverage across the mitogenome by removing reads that likely originated from numts. We compared the efficacy of Numt Parser to 2 other bioinformatic approaches that can be used to account for numt contamination. We found that Numt Parser outperformed approaches that rely only on read alignment or Basic Local Alignment Search Tool (BLAST) properties, and was effective at identifying sequences that likely originated from numts while having minimal impacts on the recovery of cymt reads. Numt Parser therefore improves the reconstruction of true mitogenomes, allowing for more accurate and robust biological inferences.

Genome report: chromosome-level draft assemblies of the snow leopard, African leopard, and tiger (Panthera uncia, Panthera pardus pardus, and Panthera tigris).

The big cats (genus Panthera) represent some of the most popular and charismatic species on the planet. Although some reference genomes are available for this clade, few are at the chromosome level, inhibiting high-resolution genomic studies. We assembled genomes from 3 members of the genus, the tiger (Panthera tigris), the snow leopard (Panthera uncia), and the African leopard (Panthera pardus pardus), at chromosome or near-chromosome level. We used a combination of short- and long-read technologies, as well as proximity ligation data from Hi-C technology, to achieve high continuity and contiguity for each individual. We hope that these genomes will aid in further evolutionary and conservation research of this iconic group of mammals.

Chromosome-level Genome Assembly of the High-altitude Leopard (Panthera pardus) Sheds Light on Its Environmental Adaptation.

The leopard (Panthera pardus) has the largest natural distribution from low- to high-altitude areas of any wild felid species, but recent studies have revealed that leopards have disappeared from large areas, probably owing to poaching, a decline of prey species, and habitat degradation. Here, we reported the chromosome-scale genome assembly of the high-altitude leopard (HL) based on nanopore sequencing and high-throughput chromatin conformation capture (Hi-C) technology. Panthera genomes revealed similar repeat composition, and there was an appreciably conserved synteny between HL and the other two Panthera genomes. Divergence time analysis based on the whole genomes revealed that the HL and the low-altitude leopard differentiate from a common ancestor ∼2.2 Ma. Through comparative genomics analyses, we found molecular genetic signatures that may reflect high-altitude adaptation of the HL. Three HL-specific missense mutations were detected in two positively selected genes, that is, ITGA7 (Ala112Gly, Asp113Val, and Gln115Pro) and NOTCH2 (Ala2398Ser), which are likely to be associated with hypoxia adaptation. The chromosome-level genome of the HL provides valuable resources for the investigation of high-altitude adaptation and protection management of the vulnerable leopard.

Efficient and cost-effective non-invasive population monitoring as a method to assess the genetic diversity of the last remaining population of Amur leopard (Panthera pardus orientalis) in the Russia Far East.

Small populations of the endangered species are more vulnerable to extinction and hence require periodic genetic monitoring to establish and revisit the conservation strategies. The Amur leopard is critically endangered with about 100 individuals in the wild. In this study, we developed a simple and cost-effective noninvasive genetic monitoring protocol for Amur leopards. Also, we investigated the impact of fecal sample's age, storage, and collection season on microsatellite genotyping success and data quality. We identified 89 leopard scats out of the 342 fecal samples collected from Land of the Leopard between 2014-2019. Microsatellite genotyping using 12 markers optimized in 3 multiplex PCR reactions reveals presence of at least 24 leopard individuals (18 males and 6 females). There was a significant difference in the success rate of genotyping depending on the time from feces deposition to collection (p = 0.014, Fisher's exact test), with better genotyping success for samples having <2 weeks of environmental exposure. Amur leopard genetic diversity was found low (Ho- 0.33, HE- 0.35, and NA- 2.57) with no visible population substructure and recent bottleneck signature. Although a historical bottleneck footprint was observed. Mitochondrial DNA diversity was also found low with two haplotypes differing by a point mutation reported in 1,769 bp of investigated sequence covering parts of cytochrome b gene (846 bp), NADH-5 gene (611 bp) and control region (312 bp). We recommend periodic genetic monitoring of wild Amur leopards following the proposed methodology to achieve cost effectiveness and efficiency.

Whole genome sequencing and the application of a SNP panel reveal primary evolutionary lineages and genomic variation in the lion (Panthera leo).

BACKGROUND: Previous phylogeographic studies of the lion (Panthera leo) have improved our insight into the distribution of genetic variation, as well as a revised taxonomy which now recognizes a northern (Panthera leo leo) and a southern (Panthera leo melanochaita) subspecies. However, existing whole range phylogeographic studies on lions either consist of very limited numbers of samples, or are focused on mitochondrial DNA and/or a limited set of microsatellites. The geographic extent of genetic lineages and their phylogenetic relationships remain uncertain, clouded by massive sampling gaps, sex-biased dispersal and incomplete lineage sorting. RESULTS: In this study we present results of low depth whole genome sequencing and subsequent variant calling in ten lions sampled throughout the geographic range, resulting in the discovery of >150,000 Single Nucleotide Polymorphisms (SNPs). Phylogenetic analyses revealed the same basal split between northern and southern populations, as well as four population clusters on a more local scale. Further, we designed a SNP panel, including 125 autosomal and 14 mitochondrial SNPs, which was tested on >200 lions from across their range. Results allow us to assign individuals to one of these four major clades (West & Central Africa, India, East Africa, or Southern Africa) and delineate these clades in more detail. CONCLUSIONS: The results presented here, particularly the validated SNP panel, have important applications, not only for studying populations on a local geographic scale, but also for tracing samples of unknown origin for forensic purposes, and for guiding conservation management of ex situ populations. Thus, these genomic resources not only contribute to our understanding of the evolutionary history of the lion, but may also play a crucial role in conservation efforts aimed at protecting the species in its full diversity.

Population genetic attributes of common leopard (Panthera pardus fusca) from Uttarkashi, Western Himalayas.

BACKGROUND: The common leopard (Panthera pardus fusca), which persists in most of its historic range, is experiencing steady population decline due to habitat loss, anthrophonic disturbances, illegal poaching for their body parts, and retaliatory killings in response to the leopard-human conflicts. METHODS AND RESULTS: We analysed 143 scats samples and identified 32 unique leopards following a selected panel of seven loci with cumulative PID sibs 5.30E-04. We observed moderate genetic diversity at nuclear (Ho = 0.600 ± 0.06) and mitochondrial markers (Hd = 0.569 ± 0.009; π = 0.001 ± 0.0002) and found sub-structuring in the leopard population at Uttarkashi, Western Himalayas. CONCLUSIONS: The present study exhibits the utility of non-invasive genetics in monitoring the leopard population and paves the path to investigate population genetic parameters in further studies.

Evidence of spatial genetic structure in a snow leopard population from Gansu, China.

Understanding the spatial structure of genetic diversity provides insights into a populations' genetic status and enables assessment of its capacity to counteract the effects of genetic drift. Such knowledge is particularly scarce for the snow leopard, a conservation flagship species of Central Asia mountains. Focusing on a snow leopard population in the Qilian mountains of Gansu Province, China, we characterised the spatial genetic patterns by incorporating spatially explicit indices of diversity and multivariate analyses, based on different inertia levels of Principal Component Analysis (PCA). We compared two datasets differing in the number of loci and individuals. We found that genetic patterns were significantly spatially structured and were characterised by a broad geographical division coupled with a fine-scale cline of differentiation. Genetic admixture was detected in two adjoining core areas characterised by higher effective population size and allelic diversity, compared to peripheral localities. The power to detect significant spatial relationships depended primarily on the number of loci, and secondarily on the number of PCA axes. Spatial patterns and indices of diversity highlighted the cryptic structure of snow leopard genetic diversity, likely driven by its ability to disperse over large distances. In combination, the species' low allelic richness and large dispersal ability result in weak genetic differentiation related to major geographical features and isolation by distance. This study illustrates how cryptic genetic patterns can be investigated and analysed at a fine spatial scale, providing insights into the spatially variable isolation effects of both geographic distance and landscape resistance.

Age estimation using methylation-sensitive high-resolution melting (MS-HRM) in both healthy felines and those with chronic kidney disease.

Age is an important ecological tool in wildlife conservation. However, it is difficult to estimate in most animals, including felines-most of whom are endangered. Here, we developed the first DNA methylation-based age-estimation technique-as an alternative to current age-estimation methods-for two feline species that share a relatively long genetic distance with each other: domestic cat (Felis catus; 79 blood samples) and an endangered Panthera, the snow leopard (Panthera uncia; 11 blood samples). We measured the methylation rates of two gene regions using methylation-sensitive high-resolution melting (MS-HRM). Domestic cat age was estimated with a mean absolute deviation (MAD) of 3.83 years. Health conditions influenced accuracy of the model. Specifically, the models built on cats with chronic kidney disease (CKD) had lower accuracy than those built on healthy cats. The snow leopard-specific model (i.e. the model that resets the model settings for snow leopards) had a better accuracy (MAD = 2.10 years) than that obtained on using the domestic cat model directly. This implies that our markers could be utilised across species, although changing the model settings when targeting different species could lead to better estimation accuracy. The snow leopard-specific model also successfully distinguished between sexually immature and mature individuals.

Whole genome survey of big cats (Genus: Panthera) identifies novel microsatellites of utility in conservation genetic study.

Big cats (Genus: Panthera) are among the most threatened mammal groups of the world, owing to hunting, habitat loss, and illegal transnational trade. Conservation genetic studies and effective curbs on poaching are important for the conservation of these charismatic apex predators. A limited number of microsatellite markers exists for Panthera species and researchers often cross-amplify domestic cat microsatellites to study these species. We conducted data mining of seven Panthera genome sequences to discover microsatellites for conservation genetic studies of four threatened big cat species. A total of 32 polymorphic microsatellite loci were identified in silico and tested with 152 big cats, and were found polymorphic in most of the tested species. We propose a set of 12 novel microsatellite markers for use in conservation genetics and wildlife forensic investigations of big cat species. Cumulatively, these markers have a high discriminatory power of one in a million for unrelated individuals and one in a thousand for siblings. Similar PCR conditions of these markers increase the prospects of achieving efficient multiplex PCR assays. This study is a pioneering attempt to synthesise genome wide microsatellite markers for big cats.

Developmental validation of Oxford Nanopore Technology MinION sequence data and the NGSpeciesID bioinformatic pipeline for forensic genetic species identification.

Species identification of non-human biological evidence through DNA nucleotide sequencing is routinely used for forensic genetic analysis to support law enforcement. The gold standard for forensic genetics is conventional Sanger sequencing; however, this is gradually being replaced by high-throughput sequencing (HTS) approaches which can generate millions of individual reads in a single experiment. HTS sequencing, which now dominates molecular biology research, has already been demonstrated for use in a number of forensic genetic analysis applications, including species identification. However, the generation of HTS data to date requires expensive equipment and is cost-effective only when large numbers of samples are analysed simultaneously. The Oxford Nanopore Technologies (ONT) MinION™ is an affordable and small footprint DNA sequencing device with the potential to quickly deliver reliable and cost effective data. However, there has been no formal validation of forensic species identification using high-throughput (deep read) sequence data from the MinION making it currently impractical for many wildlife forensic end-users. Here, we present a MinION deep read sequence data validation study for species identification. First, we tested whether the clustering-based bioinformatics pipeline NGSpeciesID can be used to generate an accurate consensus sequence for species identification. Second, we systematically evaluated the read variation distribution around the generated consensus sequences to understand what confidence we have in the accuracy of the resulting consensus sequence and to determine how to interpret individual sample results. Finally, we investigated the impact of differences between the MinION consensus and Sanger control sequences on correct species identification to understand the ability and accuracy of the MinION consensus sequence to differentiate the true species from the next most similar species. This validation study establishes that ONT MinION sequence data used in conjunction with the NGSpeciesID pipeline can produce consensus DNA sequences of sufficient accuracy for forensic genetic species identification.